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Q" 1 Abstract 



In this paper, we investigate spherically symmetric perfect fluid 
gravitational collapse in metric f(R) gravity. We take non-static 
spherically symmetric metric in the interior region and static spheri- 
^- ■ cally symmetric metric in the exterior region of a star. The junction 

conditions between interior and exterior spacetimes are derived. The 
field equations in f(R) theory are solved using the assumption of con- 
CN ' stant Ricci scalar. Inserting their solution into junction conditions, 

£ — , ■ the gravitational mass is found. Further, the apparent horizons and 

their time of formation is discussed. We conclude that the constant 
scalar curvature term f(Ro) acts as a source of repulsive force and 
thus slows down the collapse of matter. The comparison with the cor- 
responding results available in general relativity indicates that f(Ro) 
/\ ' plays the role of the cosmological constant. 
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1 Introduction 

New evidence coming from astrophysics and cosmology has revealed a quite 
unexpected picture of the expanding universe. This cosmic acceleration is 
directly supported by high red-shift supernovae, large scale structure and 
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weak lensing. What causing the repulsive effect in the cosmic expansion is 
one of the mysteries of modern era. Theoretically, the phenomenon can be 
accounted for either by supplementing the energy-momentum tensor by ex- 
otic matter component with large negative pressure (dark energy) [JJ-J3] or 
by modifying gravity itself. There exist a family of modified theories of grav- 
ity which extend in some way general relativity, for example, scalar tensor 
theory of gravity, Brans-Dick theory, string theory, Guass-Bonnet theory and 
f(R) theories of gravity. 

The f(R) theories of gravity provide the very natural gravitational alter- 
native for dark energy. It has given birth to a great number of papers in 
recent years, for example, @]-[B]. Several features including solar system test 
[7J, Newtonian limit [5], gravitational stability [9] and singularity problems 
[TO] are exhaustively discussed in this theory. General principles such as the 
so-called energy conditions have been used to place constraints on f(R) the- 
ory [11]. Amendola et al. [12] have derived the conditions under which dark 
energy f(R) models are cosmologically viable. Erickcek et al. [13] found the 
unique exterior solution for a stellar object by matching it with interior so- 
lution in the presence of matter sources in metric f(R) gravity. Kainulainen 
et al. [H] studied the interior spacetime of stars in Palatini f(R) gravity. 

Gravitational collapse is one of the most important and thorny problems 
in general relativity. It is at the heart of structure formation in the universe 
through which clusters of galaxies, stellar groups, stars and planets came into 
being. One can study the gravitational collapse by taking the interior and 
exterior regions of spacetime. The proper junction conditions are formed that 
allow the smooth matching of these regions. Oppenheimer and Synder [15] 
were the pioneers to study gravitational collapse. Herrera and Santos [T6] 
discussed the dynamics of gravitational collapse which undergoes dissipation 
in the form of heat flow and radiation. 

It has been predicted that gravitational collapse of massive objects leads 
to the formation of spacetime singularities |17j.[18|. These singularities can 
be of two types: an observable singularity (naked singularity) and a singu- 
larity that cannot be observed (black hole). Penrose [19] suggested a cosmic 
censorship conjecture (CCC) according to which the final outcome of a grav- 
itational collapse must be a black hole. Most of the studies in gravitational 
collapse have considered the Tolman-Bondi-Lemaitor spherically symmet- 
ric spacetimes containing irrotational dust to support or disprove the CCC 
[20], [21]. A general conclusion from these studies is that a central curvature 
singularity forms but its local or global visibility depends on the initial data. 



An extensive study of gravitational collapse has been carried out in the 
presence of cosmological constant. Markovic and Shapiro [22] generalized the 
pioneer work in the presence of a positive cosmological constant. Lake [23] 
extended this study with both positive and negative cosmological constants. 
Sharif and Ahmad [24J investigated gravitational collapse of a perfect fluid 
with a positive cosmological constant and generalized this work to five and 
higher dimensional spacetimes [25J . Recently, Sharif and Abbas [2H] extended 
the same work for electromagnetic charged perfect fluid. In this paper, we 
explore gravitational collapse of a perfect fluid in metric f(R) gravity. 

The paper is organized as follows. Junction conditions are discussed in 
the next section. In section 3, we solve the field equations in metric f(R) 
gravity. Here we assume the condition of constant scalar curvature which 
is very large to meet the requirement of gravitational collapse. Section 4 
is devoted to study the apparent horizons and their time of formation. We 
summarize and discuss the results in section 5. 

2 Junction Conditions 

In this section, we derive the conditions to be satisfied at the surface of a 
collapsing perfect fluid sphere. We assume spherical symmetry about an 
origin O and a 3-D hypersurface S centered at O which divides spherically 
symmetric spacetime into interior and exterior regions referred as V~ and 
V + respectively. The interior spacetime to £ can be described by the line 
element 

dsl = dt 2 - X 2 dr 2 - Y 2 (d6 2 + sin 2 6d<p 2 ), (2.1) 

where X and Y are functions of t and r. For the exterior spacetime to X, we 
take the Schwarzschild spacetime 

ds 2 + = (1 - — )dT 2 ^17 dR2 ~ R2 ( de2 + sin2 ed( ?)i ( 2 - 2 ) 

R 1_ -fT 

where M is a constant. In accordance with Israel junction conditions [28J, 
we suppose that the first and second fundamental forms inherited by S from 
the interior and exterior spacetimes are same, i.e., 

1. The continuity of line element over £ gives 

(dsi) E = (dslh = dsl. (2-3) 



2. The continuity of extrinsic curvature over £ yields 

[ICy] = K+-Kr. = 0, (j, j = 0, 2, 3) (2.4) 

where K^ is the extrinsic curvature defined as 

, ,, d 2 xl dx>idx v +s , , , 

^ = -^(^7 + r ^^f)' (<W = 0,1,2,3). (2.5) 

Here £* correspond to the coordinates on E, i^ are the coordinates in V , the 
Christoffel symbols T^ are calculated from the interior or exterior spacetimes 
and n^ are components of outward unit normals to E in the coordinates x±. 
In terms of interior and exterior spacetime coordinates respectively, the 
equation of hypersurface is given as follows 

h-(r,t) = r-r s = 0, (2.6) 

h+(R } T)=R-R< B (T) = 0, (2.7) 

where r^ is a constant. When we make use of Eq. fl2.6D in Eq. fl2.lD . the metric 
on £ becomes 

(ds 2 _)x = dt 2 -Y 2 (rx,t)(d6 2 + sm6 2 d(j) 2 ). (2.8) 

Also, Eqs.flZlD and (J2J2) yield 



#^~ 1_ 2M^^T" 



(d<)s = I 1 - ^ 5m (-FF 1 ) I dT 2 - RUd6 2 + sind'd^). (2.9) 

Here we assume that 



2A/ 



2M 1 ,dRs^ 2 



!- — -7-^M(^rl> ' ( 2 - 10 ) 

Re 



so that T is a timelike coordinate. From Eqs. fl2.3jl . ( 12. 8ft and (12.91) . it follows 

that 





i?s 


= Y(r s ,t), 


(2.11) 


2M 


1 fdRj]^ 


| ciT = dt. 


(2.12) 



The outward unit normals in regions V and V + , respectively, are found 
from Eqs.flMD and (EZD 

n- = (0,X(r E ,t),0,0), (2.13) 

n+ = (-ife,r,0,0). (2.14) 

The components of extrinsic curvature K t - become 

K 00 = 0, (2.15) 

fYY'\ 
K 22 = csc 2 ^33=^J , (2.16) 

K+ __ (fjf f B | ZM&T M(R-2M)fA 

K ™ ~ [ RT " TR+ R(R-2M) W J s ' (2 - 17) 

K+ = esc 2 ^+= [T(i2-2M)] E , (2.18) 

where dot and prime mean differentiation with respect to t and r respectively. 
From Eq. ( 12.4ft . the continuity of extrinsic curvature gives 

K+ = 0, K+ 2 = K 2 ~ 2 . (2.19) 

Using Eqs. (l2TT5|) - (l2TT9|) along with Eqs.fl2HID and d2H2J, the junction con- 
ditions turn out to be 

(XY'-XY% = 0, (2.20) 

M= (- + -Y 2 - -^Y' 2 ) . (2.21) 

\2 2 2X 2 y E v ; 

Equations fETTll . (I2TT2D . (12^201 and ( 12T211 are the necessary and sufficient 
conditions for the smooth matching of the interior and exterior regions. 



3 Field Equations in f(R) Gravity and their 
Solution 

The action of f(R) gravity in the presence of matter described by the matter 
Lagrangian Lm is given by [1] 



S — I d 4 Xy/^g 



f{R) + L 



(3.22) 



where the matter Lagrangian depends on the metric g^ u and the matter fields. 
Varying this action with respect to metric tensor, we obtain the following field 
equations 

F{R)R IW - \f{R)g^ - V„V„F(i2) + g^UF{R) = kT^, (3.23) 

where F(R) = df(R)/dR, D = V M V M with V M representing the covariant 
derivative and k(= 8n) is the coupling constant in gravitational units. These 
are the fourth order partial differential equations in the metric tensor. The 
fourth order is due to the last two terms on the left side of the equation. 
Capozziello, et al. have discussed the general solutions of these fourth order 
equation in the context of spherically symmetry [27J. If we take f(R) = R, 
these equations reduce to the field equations in general relativity. The energy- 
momentum tensor for perfect fluid is 

T^u = (p + p)u^Uu ~ V9iiv, (3.24) 

where p is the energy density, p is the pressure and w M = 5® is the four- vector 
velocity in co- moving coordinates. Using Eq. ( 13.24]) in Eq. (l3.23p . we have 

R»u = j^-[8iT((p + p)u fl u u -pg flu ) + -f(R)g llu + V ll V u F(R)-g flu nF(R)]. 

(3.25) 
Equations ( 13. 25ft for the interior spacetime take the form 

X Y 1 . 1 „ , , ,-F"(R) XF(R) 

Roo = -x- 2 Y-FjR) [8np+ 2 m - { ^^ + X 

XF'(R) , 2YF(R) 2Y'F'(R) ^ 
+ X 3 + — x^r- )], (3-26) 

X XY 2 Y" Y'X' 

R u = 2 + — ? 

X XY X 2 Y XY l 

i 1 r2,z,T>\ , 2YF (R) 2Y'F'(R). 



= ^2 f{R) - 87Tp+x{F{R) + ^r^- x*y )] > (3 - 27) 

r> / \2 XY x / \2 / \2l 

-n.22 - —y — {y) — ~xy + ^2 i~y r ^~y' ~~ ~~x ~y ~ ^y 7 

1 \ l f(m R,m + Y*(F(m F " {R) i tF{R) I kF '( R \ 

(3.28) 



-R33 
Roi 



sin 2 6R. 
V 



22- 



XV 



Y XY " F(R) 



[F> 



X 

x J 



F>] 



(3.29) 
(3.30) 



To solve this set of the field equations (I3.26p - (l3.30p . we need to integrate 
Eq. ( 13.30]) to get the explicit value of X. It follows from Eq. ()3.30p that 



X 



1Y'F + F'Y 
2Y'F + F'Y' 



dt. 



We see that this cannot be integrated unless we have some value of F such 
that derivative of the denominator should be the numerator. The condition 
of constant scalar curvature (R = R ), according to which F(R ) = constant, 
enables us to solve this integral. This means that the Ricci scalar obtained 
by contracting Eq. 03.25]) must be constant which is possible only if we take 
pressure and density constant, i.e., p = p , p = p . Consequently, the field 
equations become 



4- 2 v=Fm l& ™ +1 2 m)] ' 



X XY 

2 1 1 — 

X XY X 2l Y 



Y 
Y 



Y, 



2 Y" Y'X' 
XY 



XY 

2 Y" 



F(R ) l 2 



(3.31) 



f(Ro)-8irp }, (3.32) 



,X, 






Y' 
r h(Ro)-S7rp ], 



F(R ) l 2 
V XY' 

Y~ xY 



0. 



(3.33) 
(3.34) 



Now we solve Eqs. (13 .3 IT) - ( 13. 34jl . Integration of Eq. ( 13.341) with respect to 
t yields 

A' = Y -, (3.35) 

where W = W(r) is an arbitrary function of r. Using this equation in 
Eqs. (l33T]) -( l3T341) . it follows that 



9 k ( Y 4- - - = 

Y Y' Y 2 F{Rq) 



4tt(po - Po) - -f(Ro) 



(3.36) 



Integrating this equation with respect to t, it follows that 



m 



Y z = W z -l + 2— + 



Y F(Ro) 



4tt(po - Po) - ^f{Ro) 



Y 2 

T , (3-37) 



where m = m(r) is an arbitrary function of r and is related to the mass of 
the collapsing system. Substituting Eqs.( l3.35"j) . (J3.37P into Eq.f l3.31j) . we get 



47T 



m 



'V2 



F(R 



-(po + P o)Y'Y 



Integrating this equation with respect to r, we obtain 

, > 4:71 



m\r) 



F(Ro) 



(Po + Po) / Y'Y 2 dr + c(t), 



(3.38) 



(3.39) 



where c(t) is a function of integration. Since mass cannot be negative due to 
physical reasons, we must take the function m(r) positive. Using Eqs. ( 13.351) 
and (I3.37P into the junction condition Eq. (12.211) . it follows that 



M = m + 



QF(R ) 



4tt(po - Po) - 7>f( R °) 



Y\ 



(3.40) 



The total energy M(r, t) up to a radius r at time t inside the hypersurface 
S can be calculated by using the definition of mass function [29J given by 



M(r,t) = -Y(l+g^Y ll Y u ). 



(3.41) 



For the interior metric, it becomes 



M(r,t) = iy[l + F 2 -(^) 2 ]. 
Replacing Eqs. (13351) and (13T37J) in Eq. (ET32j) . we obtain 



M(r,t) = m(r) + 



6F(R 



4tt(po - Po) - ^f(Ro) 



Y" 



Here we assume 



F(Ro) 



Mpo - Po) - ^f(Ro) 



>0 



(3.42) 



(3.43) 



(3.44) 



and the condition 



W{r) = 1. 



(3.45) 



Using Eqs. fl3.35p . f !3.3T|) and (I3.45J) so that we have analytic solutions in 
closed form as follows 



Y 



X 



6mF(R 



4tt(po - Po) - 5/CR0) 
6mF(R ) 



sinh 3 a(r,t), 



(3.46) 



m 



A-k(pq - p ) - lf(R ) J 3m 
+ t'Jr 



[- — sinh a (r, t) 



mpo - po) - \fm cosh a(r> t)] sinh - a(r . 0) (3 47) 



3F(i? ) 



where 



, M 3(47r(p - go) - \f{R )) u 



(3.48) 



4F(i? ) 

Here t s (r) is an arbitrary function of r which denotes the time formation of 
singularity for a particular shell at distance r. 

When f(Ro) — > 87r(p —p ), the above solution corresponds to the Tolman- 
Bondi solution (30] 



lim X(r,t) 

f(Ro)-*8-rr(po-po) 



lim Y(r,t) 

f(R )->8ir(p -p ) 



m'{t s -t) + 2mt' s 


(3.49) 


[6m 2 {t s -t)}t 


9m , , 9 

~Y(ts-t) 2 


3 


(3.50) 



4 Apparent Horizons 

The apparent horizon is formed when the boundary of trapped two spheres 
is formed. The formula to find such a boundary with null outward normals 
is given as follows 

v2 r y ^2_n ( 451 ) 



g^Y^Y. = Y A - (-Y = 0. 



X 



Using Eqs. fl3T35|) and fISTBTD . we get 



F(ifc) 47T(p °-^-2 /(i?o) 



y 3 - 3F + 6m = 0. 



(4.52) 



The values of Y yield the apparent horizons. For f(Ro) = 8n(p — p ), we 
have Y = 2m, i.e., Schwarzschild horizon. Equation (I4.52p can have the fol- 
lowing positive roots. 



Case (i): For 3m < * —. v 1 jtttt , we obtain two horizons, namely 

cosmological horizon, Y c , and black hole horizon, Ybh, i.e. 



* - ^Tw -f (4 ' 53) 



* - - LTw (-f-^' (4 ' 54) 



where 



cos(p = — 3mJ — = — - — -. (4.55) 

W 47^0 - Po) - |/(i?o) l ^ 



If we take m = 0, it follows from Eqs. fl4.53p and (14.541) that 



V < = \ / T"? 3F ^\*,I,^ and ^ = °- ( 456 ) 

1, 4tt(po-Po)-|/( j Ro) 

For m 7^ and /(i?o) 7^ 87r(j?o — Po)? ^c an d ^Wi can be generalized 
respectively. 



Case (ii): For 3m = \l-r-, — °i , . , there is only one positive root which 



47r(p -po)-|/(Ro)' 

corresponds to a single horizon i.e 



n - y "-^fo-*f- */(«.)■ (4 ' 57) 

This shows that both horizons coincide. The ranges for these horizons can 
be written as follows 



U 4vr(p - Po) - 5/OR0) V 4vr(po - Po) - |/(-Ro) 

The largest proper area of the black hole horizon is given by 

^ - ^y,,^ (^9) 

Svrpobo - Po) - f(Ro) 

10 



and the cosmological horizon has its area between 



S;r/i//i!) and "^ F [ R t -, (4.60) 



87r(po-/oo)-/(i2o) " 87r(po-Po)-/(^o) 



Case (iii): For 3m > J j— ( — yZXTTr T > there are no positive roots and 



47r(p -po)-|/(flo)' 

consequently there are no apparent horizons. 

Now we calculate the time of formation for the apparent horizon using 
Eqs. (B36J and (1132]) 



'*-!,' ; ^rWr ^-l)^ (n=l,2). (4.61) 
12tt(po-Po) -|/(i2o) 2m 

When f(Ro) — >■ 87r(p — Po)> this corresponds to Tolman-Bondi solution [30] 
given as 

U = * s - -m. (4.62) 



From Eq. ( 14.61]) . we can write 



y 

— — = cosh 2 a n , (4.63) 

2m V ; 



where 



'3(47r(po-A))-|/(i2o)) 



ajr.t). { i- 4 ;^ q) 2 ^' Ur)-tn)\. (4.64) 

Equations (I4.58P and (I4.6ip imply that Y c > Y^ and t 2 > ti respectively. 
Here ti denotes the time formation of cosmological horizon and t 2 denotes 
the time formation of black hole horizon. The inequality t 2 > t\ indicates 
that the cosmological horizon forms earlier than the black hole horizon. 

The time difference between the formation of cosmological horizon and 
singularity and the formation of black hole horizon and singularity respec- 
tively can be found as follows. Using Eqs. 04.53]) - (14. 55[) . it follows that 

d(B _ 1 ( sinf + 3cosf\ <Q ^ (465) 



dm m \ simp cosip 

<& = 1 ( sm^ | 3, 
dm m \ simp cosip 

11 



rf(|k) i / sin^i 3cos^ , 

' ' + ^— > 0. (4.66) 



The time difference between the formation of singularity and apparent hori- 
zons is 

T n = t s - t n . (4.67) 

It follows from Eq.f ET63]) that 



dT n 1 



d \2m) sinha n cosha n y' F p^[47r(po - Po) - |/(^o)] 
Using Eqs.( fl~65j) and fl468|) . we get 



(4.68) 



dTi dT x d( Y < 



2m' 



dm d(^) dm m /Z[4 ir (p -p )-lf(R )]sinha 1 cosha 1 



F(Ro) 

sin f 3 cos | \ 

^ + ^ < 0(4.69) 

simp costp J 

It shows that T\ is a decreasing function of mass m. This means that time 
interval between the formation of cosmological horizon and singularity is 
decreased with the increase of mass. Similarly, from Eqs. fl4.66p and f)4.68p . 
we get 

dTo 1 



dm rnJpA-s [47r(po - Po) - |/(-Ro)] sinh a 2 cosh a 2 



F(R ) L-" U'U ruy 2 



sin (^) Scos^ 4 ^ 



+ ^— > 0. (4.70) 

sin if cos v? 

This indicates that T 2 is an increasing function of mass m showing that the 
time difference between the formation of black hole horizon and singularity 
is increased with the increase of mass. 



5 Summary and Discussion 

In this paper we have discussed different aspects of gravitational collapse in 
metric f(R) gravity. For this purpose, the field equations are solved using 
the assumption of constant scalar curvature. To meet the requirement of 
gravitational collapse, we assume that scalar curvature is very large constant 
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quantity. We obtain two physical apparent horizons (Eq. (j4.52p ) named as 
black hole horizon and cosmological horizon. It is concluded that black hole 
horizon requires more time for its formation as compared to cosmological 
horizon. However, both the horizons form earlier than singularity which 
indicates that singularity is covered, i.e, black hole. In this way, f(R) gravity 
seems to support CCC. 

Now we discuss how the constant scalar curvature term f(Ro) affects 
the rate of gravitational collapse. For the exterior metric, the Newtonian 
potential, = |(1 — goo), becomes 






MPo - Po) - 7>f(Ro) 



(5.71) 



Taking derivative of this equation with respect to R, we obtain the corre- 
sponding Newtonian force 



F = -Z + 



R 



R 2 3F(R ) 



MPo - Po) - g/ORo) 



(5.72) 



This force has no effect on collapsing matter for the following values of m 
and R 



m 



R 



^ ¥ ^ r) [Mpo-Po)- 1 2 f(Ro)} 
1 

^fe[ 47r (po-Po)-if(^)] 



However, the force will be repulsive (positive) if the values of m and R are 
larger than the above values, provided that p } R > [4-7r(p — p ) — |/(i2o)] > 0. 
In other words, the repulsive force can be generated by f(Ro) if f(Ro) < 
87r(p — Po) over the entire range of the collapsing sphere. Hence, this repulsive 
force slows down the collapse of matter. The case of attractive (negative) 
force contradicts Eq. (l3.44j) and hence our results do not favor the faster 
rate of the collapsing process. The above results can also be verified from 
Eq. (13.371) indicating the acceleration of collapsing matter as follows 



Y 2 F(R ) 



4tt(po - Po) - 7>f(Ro) 



Y 
J' 



(5.73) 
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The analysis of positive and negative acceleration provides the same results 
on f(Ro) as we obtained from the Newtonian force. 

We would like to mention here why to explore the constant curvature 
solution in a typical modified modified f(R) gravity. To generate the accel- 
erating expansion in the present universe, f(R) could be a small constant 
[22]. The universe starts from the inflation driven by the effective cosmo- 
logical constant at the early stage, where curvature is very large. When the 
curvature becomes smaller, the effective cosmological constant also becomes 
smaller. There appears the small effective cosmological constant when the 
density of the radiation and the matter becomes small and the curvature 
goes to the value Rq. Thus expansion could start and cosmological constant 
can be identified as f(Ro) in the present accelerating era. 

It is well-known that only those f(R) models are cosmo logically viable 
whose solution correspond to general relativity. The class of constant scalar 
curvature solutions with spherically symmetry contains black hole solutions 
in the presence of cosmological constant. For example, the Schwarzschild anti 
de-Sitter solutions and all the topological solutions associated with a negative 
A. e ff. Exact solutions of cylindrical symmetric spacetimes with assumptions 
of constant scalar curvature are applicable to the exterior of a string. 

Finally, we conclude that the term f(Ro) plays the same role as that 
of the cosmological constant in the Einstein field equations |24|. In general 
relativity, the cosmological constant acts as a repulsive force (hence slows 
down the collapse of matter) and same does f(Ro). It would be worthwhile 
to investigate these issues in f(R) theory by dropping the assumption of 
scalar curvature for complete understanding of gravitational collapse. 
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